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ABSTRACT 

Validation and simulations of a real-time 
dynamic cabin model were conducted on the sorbent- 
based atmosphere revitalization system for Orion. The 
dynamic cabin model, which updates the concentration 
of H 2 0 and C0 2 every second during the simulation, was 
able to predict the steady state model values for H 2 0 
and C0 2 for long periods of steady metabolic production 
for a 4-person crew. It also showed similar trends for 
the exercise periods, where there were quick changes in 
production rates. Once validated, the cabin model was 
used to determine the effects of feed flow rate, cabin 
volume and column volume. A higher feed flow rate 
reduced the cabin concentrations only slightly over the 
base case, a larger cabin volume was able to reduce the 
cabin concentrations even further, and the lower column 
volume led to much higher cabin concentrations. 
Finally, the cabin model was used to determine the 
effect of the amount of silica gel in the column. As the 
amount increased, the cabin concentration of H 2 0 
decreased, but the cabin concentration of C0 2 
increased. 

INTRODUCTION 

Recent advances [1-7] in cyclic adsorption 
process modeling have revived the interests of NASA to 
assist them in the design and development of the next 
generation sorbent-based atmosphere revitalization 
(SBAR) system for water (H 2 0) and carbon dioxide 
(C0 2 ) removal from spacecraft cabins. However, there 
are fundamental differences between Earth-bound and 
spacecraft cyclic adsorption processes. These subtle 
differences present some interesting challenges to the 
design of an SBAR system. 

In space, vacuum is readily available. So, it is 
used to lower the partial pressures of the adsorbates in 
the bed to desorb them from a nearly saturated 
adsorbent. In contrast, on Earth dilutant gases are 
cheap. So, they are generally used as purge, which also 
effectively lowers the partial pressures of the adsorbates 
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in the bed to desorb them. Although a dilutant gas can 
certainly be used in spacecraft applications, the need for 
a continuous supply of it would be exorbitantly 
expensive. As a result, the development of cyclic 
adsorption processes for commercial and spacecraft 
needs has followed different paths, with extremely large 
pressure ratios generally replacing relatively high purge 
flows to foster bed regeneration. Finally, weight on the 
spacecraft is a critical issue. To reduce weight, a 
layered bed configuration is necessary. 

The concept of a layered bed for a cyclic 
adsorption process is not new. A layered bed becomes 
of interest when components of very different 
adsorptivity need to be adsorbed from a feed stream [8]. 
The basic principle is to use a weak adsorbent at the 
inlet followed by stronger adsorbents to maximize the 
use of each adsorbent and to minimize the adsorber size 
[9]. 

One application for a pressure swing adsorption 
(PSA) process using a layered bed is air purification. 
Rege et al. [10] used a 1-bed, 4-step stripping PSA 
system to remove H 2 0 and C0 2 from an air feed. The 
adsorbents used in this study were activated alumina 
and zeolite 13X. It was shown that the H 2 0 
concentration is optimum when the bed consisted of 
75% activated alumina. The lowest C0 2 concentration 
occurred when the bed was 50% activated alumina. 
This showed that the activated alumina was needed for 
H 2 0 removal while the zeolite 13X was needed for C0 2 
removal. Reynolds et al. [1 1] showed that air purification 
can be accomplished with a 2-bed 3-step PSA system 
utilizing deep vacuum without a light reflux step. The 
maximum C0 2 removal occurred with a layering of 50% 
silica gel, 17% zeolite 13X, and 33% zeolite 5A. With 
this layering, 100% of the H 2 0 and 88.0% of the C0 2 
were removed for a TO L column. However, both of 
these air purification studies allowed the system to come 
to a periodic state with a constant feed composition. 

The new Crew Exploration Vehicle (CEV), 
Orion, will not operate at periodic conditions. Instead, 
the crew and the PSA system will experience real-time 
changes in the cabin concentrations for H 2 0 and C0 2 . 



So, why is it important to create a real-time dynamic 
cabin model? First, to gain better understanding of the 
performance of the PSA-SBAR system under typical 
dynamic conditions experienced in the cabin for the 
most diverse scenarios. Second, to foster the design, 
development, and evaluation of better PSA-SBAR 
systems. The objectives of this study are to introduce 
a real-time dynamic cabin model incorporating a PSA- 
SBAR system to simulate H 2 0 and C0 2 conditions 
inside a space craft cabin, validate the dynamic cabin 
model by showing how it agrees with a long-time steady 
state model, show simulations that display the effect of 
the feed flow rate, and cabin and column volumes on the 
removal of H 2 0 and C0 2 generated from a 4-person 
crew with a given metabolic production schedule during 
a 4-day period, and to show the effect on the removal of 
H 2 0 and C0 2 due to the amount of silica gel in the PSA 
column for a 6-person crew on a 3-day mission from 
earth to the international space station (ISS). 


PSA PROCESS DESCRIPTION 


The 2-bed 3-step stripping PSA cycle depicted in Figure 
1 includes the following cycle steps: feed (F) at just 
above atmospheric pressure (P H ), countercurrent 
depressurization (CnD), or evacuation, from P H to a 
lower pressure (P L ), and feed pressurization (FP) from 
P L to P H - While one bed undergoes the FP and F steps, 
the other bed undergoes the CnD step. The heavy 
products (H 2 0 and C0 2 ) are removed from the system 
during the CnD step, whereas the light product (mainly 
N 2 and 0 2 ) is recovered during the F step. 
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Figure 1. Schematic of the 2-bed 3-step stripping PSA cycle used in 
the parametric study with a layered bed. F = feed; CnD = 
countercurrent depressurization; FP = feed pressurization; P L = low 
pressure; P H = high pressure; Pi = intermediate pressure. 


Each bed consists of three adsorbent layers - 
silica gel, zeolite 13X, and zeolite 5A. The silica gel 
encompassed the first 50% of the bed from the feed 
end, the next 17% of the column was zeolite 13X, and 
the final 33% of the column was zeolite 5A. The silica 
gel removed the bulk amount of water, the zeolite 1 3X 
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was used to remove the remaining water and the 1 bulk 
amount of carbon dioxide, and the zeolite 5A was 
utilized for trace carbon dioxide removal. 

The periodic state process performance of this 
PSA cycle is judged by the H 2 0 and CQ 2 removal. The 
H 2 0 or C0 2 removal is defined as the amount of H 2 0 or 
C0 2 removed in the CnD step divided by the amount of 
H 2 0 or C0 2 fed to the column in the F and FP steps. 

MATHEMATICAL MODEL 


The assumptions and equations used by 
Reynolds et al. [11] in their rigorous multilayer PSA 
simulator are the same utilized in this study. However, 
some modifications were necessary to accommodate 
the extension to a cabin model. A schematic of the 
cabin model is shown in Figure 2. After each second the 
model calculates new amounts, in moles, of H 2 0 and 
C0 2 present in the cabin assuming perfect mixing of the 
moles coming from the PSA column into the cabin. 

^i,j ^i,j - 1 ^i, effluent ^ astronaut ' ^ i, astronaut (^) 


where «. . is the new molar amount of component i, 
n i M is the previous molar amount of component i, 
n t fed is the molar amount of component i fed to the PSA 
system since the last update, n ieffluent is the molar 

amount of component i returned to the cabin in the 
effluent from the PSA system since the last update, 

N astronaut is the number of astronauts on-board the 

Orion, and n i astronaut is the molar amount of component i 

produced by each astronaut based on the metabolic 
production of H 2 0 and C0 2 during the activity that the 
astronaut is currently undergoing. These molar amounts 
are then converted to partial pressures using the ideal 
gas law. The partial pressure of oxygen (0 2 ) is held 
constant at its atmospheric partial pressure, while the 
amount of nitrogen (N 2 ) added to the cabin is adjusted to 
ensure that the cabin pressure remains at tho set 
pressure. Then, the molar fraction is calculated by 


n TJ 


The total number of moles (n Tj ) is calculated by 

summing the number of moles of all four components 
(H 2 0, C0 2 , 0 2 and N 2 ). Finally, these new 

concentrations are used as the molar fractions of the 
feed to the column. 

RESULTS AND DISCUSSION 

The cabin model was first validated against a 
steady state model that was based only on the total flow, 
the fraction of time that one of the columns was 
undergoing the F step, the metabolic rate of the 




astronauts, and the removal efficiency of the column. 
Figure 3 shows typical metabolic production rates during 
periods of rest, normal activity, exercise and normal 
activity for H 2 0 and C0 2 over a 24 hr period. The 
schedule is repeated for each day during a mission. As 
can be seen in the figure, there are two distinct regions 
of exercise. The first region is a period of aerobic 
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Figure 2. Schematic of the real-time dynamic cabin model used in the 
PSA-SBAR air purification study. 

exercise, while the second region is a period of strength 
training. The PSA-SBAR system was in the feed step 
for 420 seconds out of a half cycle of 480 seconds, thus 
the fraction fed was 0.875. For 10 L columns in the 
PSA-SBAR system, the dynamic cabin model had 
average removal efficiencies of 100% for H 2 0 
throughout the entire run and 83.0%, 82.9%, 79.6%, 
77.6% and 79.4% for C0 2 during periods of rest, normal 
activity, aerobic exercise, strength exercise and normal 
activity, respectively. As shown in Figure 4, the cabin 
model was able to predict the steady state model values 
for H 2 0 and C0 2 for long periods of steady metabolic 
production. It also showed similar trends for the 
exercise periods, where there were quick changes in 
production rates. 

Since the cabin model was validated against 
steady state predictions, a study was conducted to 
display the effect of the feed flow rate, and cabin and 
column volumes on the removal of H 2 0 and C0 2 
generated from a 4-person crew with a given metabolic 
production schedule for a 4-day mission. The base case 
bed characteristics and operating parameters are shown 
in Table 1. The high pressure is set at atmospheric 
pressure, and the system is run with a high-to-low 
pressure ratio of 1,000. The crew size, initial feed 
composition, feed flow rate, feed and ambient 

temperatures, and overall heat transfer coefficient were 
provided by NASA. The cycle step times are similar to 
those used in a previous study [11], and the mass 
transfer coefficient for H 2 0 on silica gel was 0.002 s' 1 . 

The effects of flow rate (Q), cabin volume 

(Vcab), and column volume (Vcol) on cabin H 2 0 and 
C0 2 concentrations (Figure 5) and PSA removal 

efficiencies (Figure 6) are illustrated with a 4-day 

mission for a 4-person crew with a H 2 0 and C0 2 


metabolic production schedule shown in Figure 3. The 
bold and thin lines correspond to H 2 0 and C0 2 , 
respectively. The dashed lines indicate the desired 
bounds (6 and 17.5 torr) for the H 2 0 cabin 
concentration. These bounds are set to eliminate any 
harm to the astronauts. If the H 2 0 
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Figure 3. Typical metabolic production rates during periods of rest, 
normal activity, exercise and normal activity over a 24 hr period for a 4- 
person crew. Line styles: bold - H 2 0; thin - C0 2 . 


Table 1. Bed characteristics and operating parameters utilized in the 
4-person crew PSA-SBAR air purification study. 


Bed Length(L) 
Bed Radius (rt>) 
Cabin Volume (Vcab) 
High Pressure (P H ) 
Low Pressure (Pl) 
Crew Size 

Initial Feed Mole Fraction: H 2 0, C0 2 , 

0 2 , N 2 

Feed Flow Rate (Qf) 
Feed Temperature (Tf) 
Ambient Temperature (T 0 ) 
Overall Heat Transfer Coefficient (h) 
Step Times (t s ): F,CnD, FP 


0.244 m 
0.081 m 
11 m 3 

101.325 kPa 
0.101325 kPa 
4 

0.0125, 0.0021 
0.2083, 0.7771 
700.00 SLPM 
291.48 K 
291.48 K 

5.68x1 O' 4 kJ/m 2 s K 
420, 480, 60 s 



Figure 4. Validation run for the dynamic cabin model utilizing 10 L 
columns for the PSA-SBAR system. 
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Figure 5. H 2 0 and C0 2 cabin concentrations from simulations of the dynamic 
cabin model for a) base case, b) effect of flow rate (Q F = 1000 SLPM), c) effect 
of cabin volume (V ca b = 44 m 3 ) and d) effect of column volume (V co i = 2.5 L). 
Line styles: bold - H 2 0; thin - C0 2 ; dashed - H 2 0 bounds; dotted - C0 2 upper 
bound. 
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Figure 6. H 2 0 and C0 2 removal efficiencies from simulations of the 

dynamic cabin model for a) base case, b) effect of flow rate (Q F = 1000 
SLPM), c) effect of cabin volume (V ca b = 44 m 3 ) and d) effect of column 
volume (Vcoi = 2.5 L). Symbols: squares - H 2 0; triangles - C0 2 . 






cohcentration goes below the lower bound, the 
astronauts’ skin will begin to crack. If the upper bound is 
exceeded, H 2 0 will begin to condense in the cabin. The 
dotted line indicates the safe upper bound (7.6 torr) for 
the C0 2 cabin concentration. If this bound is exceeded 
for an extended amount of time, it would be harmful to 
the astronauts. 

Figure 5(a) shows the base case condition. It 
can be seen that during periods of rest and normal 
activity there was too little H 2 0. At the same time the 
C0 2 was well below its bound. During these periods, 
the column was removing too much water and drying the 
cabin out. This can be avoided by slowing the feed flow 
rate down. However, during exercise periods the upper 
bounds for both H 2 0 and C0 2 were exceeded. As 
shown in Figure 6(a), the removal efficiency for both H 2 0 
and C0 2 fell shortly after the exercise period began. 
The H 2 0 removal efficiency falling meant that the H 2 0 
penetrated through the zeolite layers and returned to the 
cabin. Since the C0 2 was removed only by the zeolite 
layers, its removal efficiency fell as the H 2 0 was 
preferentially adsorbed by the zeolites. However, an 
interesting observation of Figure 5(a) was that the H 2 0 
and C0 2 concentrations returned to the pre-exercise 
level within a few hours. 

Figure 5(b) shows the H 2 0 and C0 2 
concentrations for the effect of increasing Q to 1000 
SLPM. Higher flow rates should lead to lower H 2 0 and 
C0 2 cabin concentrations. The concentrations of H 2 0 
and C0 2 in the cabin were only slightly lower than those 
from the base case. A large decrease in the cabin 
concentrations was not observed due to the lower 
removal efficiencies of the PSA unit as shown in Figure 
6(b). 

Figure 5(c) shows the effect of increasing the 
Vcab to 44 m 3 . The H 2 0 and C0 2 cabin concentrations 
were reduced, even though the removal efficiencies 
were similar to the base case. The concentrations, 
however, still exceeded the upper bounds. But, this did 
not happen until much later in the exercise; period. The 
larger cabin led to a slower response, which dampened 
the large variations observed in Figure 5(a). 


Table 2. Bed characteristics and operating parameters utilized in the 
6-person crew PSA-SBAR air purification study. 


Bed Length (L) 

0.244 m 

Bed Radius (r b ) 

0.081 m 

High Pressure (P H ) 

101.325 kPa 

Low Pressure (P L ) 

0.101325 kPa 

Crew Size 

6 

Initial Feed Mole Fraction: H 2 0, C0 2 j 

0.0125, 0.0033, 

0 2 , N 2 

0.2083, 0.7759 

Feed Flow Rate (Q F ) 

849.51 SLPM 

Feed Temperature (T F ) 

294.26 K 

Ambient Temperature (T 0 ) 

294.26 K 

Overall Heat Transfer Coefficient (h) 

5.68x1 O' 4 kJ/m 2 s K 

Step Times (t s ): F,CnD, FP 

420, 480, 60 S 


Figure 5(d) shows the effect of decreasing V C0 | 
to 2.5 L. The cabin concentrations were dramatically 


higher than that observed for the base case. In fact, the 
C0 2 concentration never came to a periodic state. The 
smaller column led to lower removal efficiencies 
resulting in higher H 2 0 and C0 2 cabin concentrations. 

The second study performed on this system 
consisted of a 6-person crew undertaking ground 
operations and a 3-day mission to dock the Orion with 
the ISS. Ground operations consisted of a 6 hour period 
of time prior to launch. Throughout the ground 
operations the PSA-SBAR system was run at 
atmospheric conditions with one system undergoing the 
F step and the other system undergoing a purge step 
that consisted of N 2 flowing in a countercurrent direction 
to remove H 2 0 and C0 2 from the adsorbent layers. 
After ground operations, there was essentially no H 2 Q or 
C0 2 on the adsorbents due to the N 2 purge. 

Table 2 shows the conditions used for the 5 
simulations run in this study. The initial feed mole 
fraction is the concentration of the components in the 
cabin at the start of ground operations. The mass 
transfer coefficient for H 2 0 on silica gel was 0.029 s' 1 . 
This number was used as a conservative estimate of 
results from recent experiments not shown here. The 
metabolic production schedule is shown for one day in 
Figure 7. Again, this schedule is repeated every day 
during the mission. Since there are 6 crew members 
and a short mission, no exercise is required. 
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Figure 7. Typical metabolic production rates during periods of normal 
activity, rest and normal activity over a 24 hr period for a 6-person 
crew. Line styles: bold - H 2 0; thin - C0 2 . 

Figure 8 shows the H 2 0 and C0 2 cabin 
concentrations for ground operations followed by a 3- 
day mission to the ISS, while Figure 9 shows the 
removal efficiencies of H 2 0 and C0 2 and the loss of 0 2 
and N 2 to space for only the 3-day mission. The first 6 
hours of Figure 8 correspond to the ground operations. 
It can be seen that the H 2 0 and C0 2 came to a periodic 
state in the cabin and that the levels in the cabin were 
well within the bounds for each component. The H 2 0 
concentration in the cabin fell below the lower bound 
during the rest period for each of the 3 days. This was 
due to the high flow rate being held constant throughout 
the entire run. To keep the H 2 0 level above the lower 
bound, the flow rate can be reduced during periods of 



rest. However, the C0 2 concentration went above its 
safe upper limit during the second day of normal activity. 
This was dqe to H 2 0 penetrating into the zeolite layers. 
Water preferentially adsorbed to the zeolite, thus C0 2 
was not be able to adsorb and be removed from the feed 
stream. This can be seen in Figure 9 as the C0 2 
removal efficiency continuously declined until about 40 
hours into the 3-day mission. One last observation from 
Figure 8 is that the H 2 0 concentration began to increase 
in the cabin after approximately 58 hours of operation 
(52 hours into the mission). This was due to H 2 0 
breaking through the column and returning to the cabin 
in the PSA effluent stream. This observation 
corresponded with the drop in the removal efficiency of 
H 2 0 shown in Figure 9. 
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Figure 8. H 2 0 and C0 2 cabin concentrations for ground operations 
and a 3-day mission to the ISS for a system containing 50% silica gel. 
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Figure 9. H 2 0 and C0 2 removal efficiencies and 0 2 and N 2 losses for 
a 3-day mission to the ISS for a system containing 50% silica gel. 
Symbols: squares - H 2 0 removal efficiency; triangles - C0 2 removal 
efficiency; empty squares - 0 2 loss; empty triangles - N 2 loss. 


The effect of the amount of silica gel in the PSA- 
SBAR column on the H 2 0 and C0 2 cabin concentrations 
were difficult to see in a graph similar to Figure 8. There 
were only slight differences between the cabin 
concentrations achieved with 35% and 55% silica gel. 
However, the differences can be seen with graphs 
similar to Figure 9. Thus, Figures 10 and 11 show the 
effect of the amount of silica gel in the PSA-SBAR 
column on the removal efficiencies of H 2 0 and C0 2 . 

As the silica gel percentage increased in the 


I 

column, the H 2 0 removal efficiency remained higher 
throughout the run as shown in Figure 10. Since the 
silica gel was used as the bulk H 2 0 remover, more silica 
gel in the column led to more H 2 0 being removed. This, 
in turn, meant that less H 2 0 needed to be adsorbed by 
the zeolite layers. Thus, breakthrough of the H 2 0 
occurred later in the run for the columns containing more 
silica gel. Therefore, the H 2 0 concentration in the cabin 
was slightly lower for the columns containing more silica 
gel. 



Figure 10. H 2 0 removal efficiency for a 3-day mission for the CEV to 
dock with the ISS. Symbols: squares - 35% silica gel; triangles - 40% 
silica gel; circles - 45% silica gel; empty squares - 50% silica gel; 
empty triangles - 55% silica gel. 

Figure 11 shows the C0 2 removal efficiency 
throughout the 3-day mission as a function of the silica 
gel percentage. The silica gel had an opposite effect on 
the C0 2 as it did on the H 2 0. For C0 2 , the removal 
efficiency decreased as the silica gel percentage 
increased. This was caused by the fact that as the silica 
gel percentage increased the percentage of zeolite 13X 
decreased. Since the C0 2 was removed only by the 
zeolite layers, the capacity to remove the C0 2 from the 
feed stream was decreased as the zeolite 13X amount 
was decreased. Therefore, the cabin concentration of 
C0 2 increased as the silica gel percentage increased. 



Figure 11. C0 2 removal efficiency for a 3-day mission for the CEV to 
dock with the ISS. Symbols: squares - 35% silica gel; triangles - 40% 
silica gel; circles - 45% silica gel; empty squares - 50% silica gel; 
empty triangles - 55% silica gel. 










CONCLUSIONS 

A rigorous cabin model was developed to aid in 
the design and development for the SBAR system on 
the new CEV, Orion. This cabin model was able to 
update the cabin concentrations of H 2 0 and C0 2 every 
second throughout the run. It was validated against a 
steady state model that was based only on the total flow, 
the fraction of time that one of the columns was 
undergoing the F step, the metabolic rate of the 
astronauts, and the removal efficiency of the column. 
The dynamic cabin model was able to predict the steady 
state model values for H 2 0 and C0 2 for long periods of 
steady metabolic production. It also shows similar 
trends for the exercise periods, where there are quick 
changes in production rates. 

Once validated, the dynamic cabin model was 
used to illustrate the effect of flow rate, cabin volume 
and column volume on the cabin concentrations of H 2 0 
and C0 2 for a 4-person crew over a 4-day period. 
Increasing the flow rate led to slightly lower cabin 
concentrations of H 2 0 and C0 2 since more of these 
components are able to interact with the adsorbent. The 
expected decrease in the cabin concentrations was not 
observed as the removal efficiency for both H 2 0 and 
C0 2 was much lower than the base case. 

An increase in the cabin volume caused a 
reduction in the H 2 0 and C0 2 cabin concentrations, 
even though the removal efficiencies were similar to the 
base case. The larger cabin led to a slower response, 
which dampened the large variations observed in the 
base case. A smaller column adsorbed less H 2 0 and 
C0 2 than the base case which resulted in dramatically 
larger cabin concentration than the base case. In fact, 
the C0 2 concentration never came to a periodic state. 
The smaller column led to lower removal efficiencies 
resulting in higher H 2 0 and C0 2 cabin concentrations. 

Finally, the cabin model was used to observe 
the effect of the silica gel percentage in the column on 
the concentrations of H 2 0 and C0 2 . As the silica gel 
percentage increased in the column, the H 2 0 removal 
efficiency remained higher throughout the run. Thus, 
breakthrough of the H 2 0 occurred later in the run for the 
columns containing more silica gel. Therefore, the H 2 0 
concentration in the cabin was slightly lower for the 
columns containing more silica gel. The silica gel, 
though, had an opposite effect on the C0 2 as it did on 
the H 2 0. For C0 2 , the removal efficiency decreased as 
the silica gel percentage increased. This was caused by 
the fact that as the silica gel percentage increased the 
percentage of zeolite 13X was decreased. Thus, the 
capacity to remove the C0 2 from the feed stream was 
decreased as the zeolite 13X amount was decreased. 
Therefore, the cabin concentration of C0 2 increased as 
the silica gel percentage increased. 
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